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Abstract 



According to previous work on magnetic monopoles, static regular solu- 
tions are nonexistent if the vacuum expectation value of the Higgs field rj is 
larger than a critical value r/cr; which is of the order of the Planck mass. In 
order to understand the properties of monopoles for r] > r/cr, we investigate 
their dynamics numerically. If r/ is large enough ?/cr)) a monopole expands 
exponentially and a wormhole structure appears around it, regardless of cou- 
pling constants and initial configuration. If is around r/cr, there are three 
types of solutions, depending on coupling constants and initial configuration: 
a monopole either expands as stated above, collapses into a black hole, or 
comes to take a stable configuration. 
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I. INTRODUCTION 



In recent years static and spherically symmetric solutions of the Einstein- Yang-Mills- 
Higgs system have been intensively studied in the literature [0-0]. One purpose of such 
investigation has been to understand the nature of black holes, especially in the context 
of the no-hair conjecture; it was shown that non-trivial black holes are stable and hence 
the monopole black hole could be one of the most plausible counterexamples. The other 
interest has been in the properties of particle-like solutions; it was shown that such regular 
monopoles exist only if the vacuum expectation value of the Higgs field rj is less than a 
critical value rj^r, which is of the order of the Planck mass mpi. This result naturally gives 
rise to the next question: what is the fate of monopoles for t] > 

Because the only static solution for rj > rjcr is the Reissner-Nordstrom black hole, we 
can expect that a monopole which is regular initially evolves into the Reissner-Nordstrom 
black hole. Even if this speculation is reasonable, we still do not know how the black hole 
formation occurs. One could imagine two alternatives: a monopole just shrinks, or its core 
continues to expand inside the black-hole horizon, just as a "child universe" [Q. 

Linde and Vilenkin independently pointed out the latter possibility in the context of 
the "topological inflation" model 0. They claimed monopoles as well as other topological 
defects expand exponentially if 77 > O(mpi). Their discussions for the Einstein-Higgs system 
were verified by our numerical simulation in |^: we found that domain walls and global 
monopoles inflate if and only if ~ 0.33mpi. The next question on this monopole inflation 
is similarly what happens to magnetic monopoles in the Einstein- Yang-Mills-Higgs system. 
Because we cannot find an answer to the question only by analyzing static solutions, we 
investigate dynamic monopole solutions in this paper. 

The plan of this paper is as follows. In Sec. II, we derive the basic equations and 
explain how we solve those dynamical equations numerically. In Sec. Ill, we offer analytic 
discussions and numerical results. Sec. IV is devoted to summary and discussions. In this 
paper we use the units c = h = 1. 



II. BASIC EQUATIONS 



The 5*0(3) Einstein- Yang-Mills-Higgs system is described by 




(2.1) 



with 
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where A'^ and are the SU{2) Yang-Mills field potential and its field strength, respectively. 

is the real triplet Higgs field and V{^) is its potential. A and e are the Higgs self coupling 
constant and the gauge coupling constant, respectively. and are the spacetime 
covariant derivative and the totally covariant derivative, respectively. The variation of 
with respect to g^^, and A" yield the Einstein equations. 



and the equations for the matter fields 
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We assume a spherically symmetric spacetime and adopt the coordinate system, 

ds" = -dt" + A^{t, r)dr^ + B^{t, ry{de^ + sin^ Odip^). (2.7) 

For the matter fields, we adapt the 't Hooft-Polyakov ansatz in such a way that we can 
apply it to a time- dependent curved spacetime: 



= ^{t, = (sin 9 cos Lp, sin 9 sin Lp, cos 9). 
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With the metric ( p. 71) and the 't Hooft-Polyakov ansatz (|2.8|) and (|2.9|) , we write down the 
field equations (|2^ )- (|2.6| ) as 
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where the overdot and the prime denote the partial derivative with respect to t and r, re- 
spectively. We have introduced the extrinsic curvature tensor of a t = constant hypersurface, 
Kij, whose components are given by 
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and we have denoted its trace hy K = Kl. 

As an initial configuration of the matter fields, we adopt the functional form of the static 
solution in a fiat spacetime with A = 0: 



$(t = 0,r) = <l>fl,t — 

W{t = 0,r) = Wflat( — 
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where c$ and are the initial size parameters. The configurations of $flat('^) and of Wfiat(^) 
are illustrated in Fig. 1. As to the time-derivative, we suppose $(t = 0, r) = w{t = 0, r) = 0. 



In order to set up consistent initial data, we have to solve the constraint equations (p.lO|) 
and ( p.ll|) . At this point, there are four unknown variables. A, B, K and K^, in the two 
constraint equations; two of the variables are arbitrarily chosen. One of the methods which is 
usually adopted is to assume K = const and A = B. In this system, however, the condition 
of i^' = const 7^ is not appropriate because the far region is asymptotically fiat. Further, in 
the range where there exists no static solution, we cannot fix K = even momentarily. As 
an alternative, thereby, we suppose A{t = 0,r) = B{t = 0,r) = 1 and solve the constraint 
equations ( ^.l(j| ) and ( |2.11| ) to determine K{t = 0,r) and K^t = 0,r). This treatment is 
suitable for this system because we obtain 
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which approaches zero as r increases; we can construct an asymptotically fiat spacetime 
without iterative integration. We have also assumed K{t = 0,r) < 0: every point in 
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the spacetime is locally expanding. The numerical boundary is fixed at r = 30 /{erj) ^ or 
60/(er/)-i. 

In order to solve the dynamical equations, we use a finite difference method with 2000 to 
10,000 meshes. Now we have six dynamical variables: A, B, K, $ and w. Equations 
( |^.12| )-( PTT3| ) provide the next time-step of A, B, K, $ and w, respectively. At each step, 
we integrate ( |2.11| ) in the r-direction to obtain In this way we have reduced spatial 
derivatives appearing in the equations, which may become seeds for numerical instability. 
The Hamiltonian constraint equation ( ^.1U| ) remains unsolved during the evolution and is 
used for checking numerical accuracy. We stop numerical computation when some errors 
exceed a few percent. 

In order to understand the spacetime structure from the numerical data, it is useful to 
observe the signs of the expansion of a null geodesic congruence. Nambu and Siino [0] also 
utilized this tool to study wormhole formation in a singlet scalar field system. For the metric 



(p.7|), the expansion 9± is written as 



0-1- — + ^±;3 



-A1±^J, (2.18) 

where k'^ = (—1, ±A~^,0,0) is an outgoing (+) or ingoing (— ) null vector. We observe the 
signs of G-t at all points in the numerical spacetime. For later convenience, we define "RI" 
as the region where both and 6~ are positive, and "RII" as the region where they are 
negative. We can interpret that the region around RI is de Sitter-like and the region around 
RII is Schwarzschild-like. The two-surface which bounds RI or RII is called an "apparent 
horizon". Later, we will use the term "black hole horizon" to refer to any boundary of 
RII. And, we will use the term "cosmological horizon" in the sense that no information 
beyond it reaches the center of a monopole; only the innermost boundary of RI is called the 
cosmological horizon. 



III. NUMERICAL RESULTS 

Before we move on to numerical simulation, we offer a rough discussion on the effect of 
the gauge fields on the gravitational field. In our previous paper 0, we investigated the 
effect for static monopole solutions and found that the gauge fields generate an attractive 
force. This property is also described in the time-dependent coordinate system ( p.7| ) as 
follows. 

In a homogeneous and isotropic spacetime, the evolution of the scale factor a{t) follows 

"-(p + ip), (3.1) 
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where p and p are the energy density and the pressure of a matter, respectively. Equation 
indicates that the sign of p + 3p determines whether the acceleration of the cosmic 
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expansion is positive or negative. We can extend this discussion to general spacetimes: 
the sign of p + Spj = —T^ + T? determines whether a local region expands with positive 



acceleration or not. A corresponding equation in the present system to ( p.lD is ( |2.12| ). At 
the origin, ( |2.12D reduces to 
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In the case of global monopoles, the second term of the right-hand side disappears, and 
hence the central region is always locally de Sitter spacetime. If gauge fields exist, however, 
the local acceleration at the center also depends on the second term. Although the exact 
value of d'^w / {AdrY cannot be determined without solving the full dynamical equations, we 
can estimate its order by use of the static solution in a flat spacetime with A = 0. Assuming 
d'^w / {AdrY = b w'^^^^{r) with b = 0(1), we have 

P + Sp.Uo = ^(-^ + ^&), (3.3) 

which lets us understand how the local expansion of the spacetime in the center depends 
on A/e^. We see that, if A/e^ -C 1, the monopole core is an attractive spacetime; while, if 
A/e^ ^ 1, it is repulsive like de Sitter spacetime. Of course, if the initial conflguration is 
quite different from that of the static solutions, i.e., b ^ 0(1), the above discussion is not 
true. The dynamics may also depend on initial conflguration. 

In what follows, by use of the method in Sec. II, we will numerically integrate the 
fleld equations (|2.1CI|) - (p?T^) . To show the results, we deflne X as a proper distance in the 

radial direction: X = Adr. We also deflne the boundaries of a monopole in two ways: 

Jo 

Xij,(t) = X at the position of $ = ri/2 and X^(t) = X at the position oi w = 1/2. We 
normalize time and length by the horizon scale deflned as Hq^ = (8nV{0) /3mpi)^^ . 

First, we check our numerical code by solving the equations for the case of weak gravity. 
In Fig. 2 we set t] = O.lmpi and A/e^ = 0.1, and give two initial conflgurations: c$ = = 1 
and 0.5. We plot the trajectories of X$(t). We flnd that the flelds behaves stably; this 
reasonable result indicates that our numerical code works well. 

From now on, we concentrate on the parameter range where no static solution exists. In 
Fig. 3 we set r] = 0.4mpi and A/e^ = 0.1, and give two initial conflgurations: c$ = = 1 
in (a) and c$ = = 10 in (b). In Fig. 3(a)(b) we plot the trajectories of X^(t) and Xw{t) 
as well as apparent horizons. The dynamics in these two cases contrast sharply: in (a) a 
monopole shrinks and black-hole horizons appear, while in (b) a cosmological horizon exists 
from the beginning and a monopole continues to expand. We also draw the distributions of 
p+ T,pi in (c) and in (d), which correspond to the results in (a) and in (b), respectively. In 
(c) the values around the center become negative at the beginning, but they bounce back to 
positive values, which conflrms that the monopole core never inflates. On the other hand. 
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in (d) the values of p + Epj around the center remain negative from the beginning. This 
behavior indicates that exponential expansion really occurs inside the monopole. These two 
results tell us that monopoles for r] > rjcr tend to be dynamic, and their dynamics depends 
on the initial configuration, contrary to the case of global monopoles. 

As we will see soon, for larger rj, monopoles are more likely to inflate rather than shrink. 
We show an example for larger rj in Fig. 4; we set rj = 0.55mpi, A/e^ = 0.1 and c$ = c^i, = 1. 
In Fig. 4(a) wc plot the trajectories of X$(t) and Xyu{t) as well as apparent horizons. 
(Please also refer to Fig. 5, which is a schematic sketch of the spacetime structure.) From 
the beginning there arc two apparent horizons, 5*1 and 5*2: 5*1 is the cosmological horizon. 
Later other two apparent horizons, S3 and 54, appear, and then S2 and 54 approach each 
other. These surfaces turn out to be black-hole horizons, S2' and 54', the moment they 
intersect. In Fig. 4(b) we draw the distributions of p + Spj. Contrary to the case of the 
contracting monopole in Fig. 3(a), the values around the center are initially positive, but 
they become negative. This suggests that, if rj is large enough, a monopole begins to expand 
exponentially even if its initial size is not so large. We also show in Fig. 4(c) the relation 
between the proper distance along the radial direction and the circumference radius, which 
indicates a wormhole structure really appears. Figure 4 lets us understand how the wormhole 
is created. Because the expanding core is causally disconnected from the outer region, such 
an isolated region is called a "child universe" . 

One may think that if > rjcr, a monopole either expands or collapses, as shown in 
Fig. 3 or 4. However, we find some cases where a monopole neither expands nor collapses. 
An example of such solutions is shown in Fig. 6 {rj = 0.3mpi and A/e^ = 1). Setting 
c$ = c^t, = 1, we show the evolution of $ in (a) and that of w in (b), and the trajectories 
of X<s>{t) and of X^(t) in (c). Although some oscillations remain outside the monopole, the 
core of the monopole approaches a stable configuration. Wc change the initial size in Fig. 
6(d), finding monopoles with any initial size behave stably. These results indicate that there 
exist stationary solutions. 

In order to see if such stable monopoles are really created in an expanding universe, 
we consider a different type of initial configurations: to give a small perturbation on the 
symmetric state (de Sitter spacetime). Specifically, we adopt a form. 



where we fix d = 0.1 and c = 20 in our analysis. In this way we see how a monopole 
configuration is formed from the nearly symmetric state. We assume A/e^ = 10 and 77/mpi = 
0.3, as is the case is Fig. 6. Figure 7 shows that, once a monopole configuration is formed, it 
approaches a stable configuration instead of continuing to expand. This result supports the 
existence of stationary solutions. The existence of such solutions looks surprising, because 
one may expect that all solutions in the parameter range where no static solution exists 
must be dynamical. The existence of "stationary" solutions, however, does not contradict 




w{t = 0,r) = 1, 



(3.4) 
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the nonexistence of "static" solutions: the solutions in Figs. 6 and 7 cannot be described with 
a static coordinate system because the size of the monopole is greater than the cosmological 
horizon. 

Finally we systematically survey the dynamics of monopoles for 0.05 < rj/rripi < 0.55 
and 0.1 < A/e^ < 10 and summarize the solutions in the A/e^-r^/mpi plane of Fig. 7. A 
square denotes a stable solution, as is the case in Fig. 2 or 6. A cross denotes the case 
where a monopole shrinks, as is the case in Fig. 3(a). A circle denotes the case where a 
monopole inflates and the wormhole structure appears, as is the case in Fig. 4. A dotted 
line indicates the maximum values of rj/m-p\ versus A/e^, depicted approximately by use of 
Fig. 6 in [Q. We vary c$ and c^, from 1 to 10, and hence some parameter points are labeled 
as two symbols. We interpret these results as follows. In the case of A/e^ > 1, a monopole 
expands exponentially if ~ 0.35mp;; this critical value has little dependence on A/e^ and 
initial configuration, and almost agrees with that for global monopoles 0. This agreement 
is quite reasonable because the effect of the gauge fields is smaller as A/e^ is larger. Below 
the critical value, a monopole tends to take a stable configuration even in the theories where 
static solutions are nonexistent. In the case of A/e^ < 1, the dynamics also depend on A/e^ 
and initial configuration. In some cases the effect of the gauge fields becomes dominant and 
a monopole shrinks and becomes a black hole. These results are consistent with our analytic 
discussions at the beginning of this section. 



IV. SUMMARY AND DISCUSSIONS 

We have studied the dynamics of magnetic monopoles numerically. Our main purpose 
has been to understand the behavior of monopoles in the case where static solutions are 
nonexistent. 

If 7] is large enough ^cv)i a monopole inflates and a wormhole structure appears around 
it. We have shown how the wormhole connected with a child universe is created. We should 
emphasize that a child universe can be generated without fine-tuned initial conditions in this 
model, contrary to the case of a trapped false vacuum bubble 0]. In the case of A/e^ > 1, 
the condition of inflation is ~ 0.35mp;, which has little dependence on A/e^ and initial 
configuration. Below the critical value, a monopole tends to take a stable configuration 
even in the theories where static solutions are nonexistent. This is true for any initial 
configuration, which indicates the existence of stationary solutions. While, in the case of 
A/e^ < 1, the dynamics also depends on A/e^ and initial configuration. In some cases the 
effect of the gauge fields becomes dominant and a monopole collapses into a black hole. 

We should notice that the condition of inflation was also estimated analytically in a 
simplified model by Tachizawa et al. They discussed the global structure of a space- 
time by regarding the inside the monopole core as de Sitter spacetime and the outside as 
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Reissner-Nordstrom spacetime. They showed that the surface of the monopole core exceeds 
a cosmological horizon if > mpi/-\/37r ^ 0.33mpi. This condition almost agrees with the 
condition of inflation for most cases in our analysis; this agreement suggests that our nu- 
merical results are reasonable as well as that their simplified model is a good approximation 
in most cases. When the effect of the gauge fields is dominant to that of the Higgs field, 
however, the spacetime is not de Sitter-like, and then the validity of the simplified model is 
lost. 

Our results as a whole support the discussions of Linde and Vilenkin who pointed 
out the possibility of monopole inflation. Actually, we have found that inflation happens in 
most cases of t] > rj^r- Further, if the initial size of a monopole is large enough, the effect of 
the gauge fields is not important, as Linde mentioned. What we have clarified more about 
this subject is there are some cases where static solutions are nonexistent but monopoles do 
not continue to expand, as the results in Figs. 6 and 7. Although we did not state this fact 
in our previous paper P] , it is also true for global monopoles. 
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Figure Captions 



FIG.l. Configurations of $flat('^) and of wnatir) in ( ^2.161 ), which are adopted as initial 
conditions for time-evolutions. 

FIG. 2. Dynamics of a monopole for the case of weak gravity. We set r] = O.lmpi and 
A/e^ = 0.1, and give two initial configurations: c$ = = 1 and 0.5. We plot the trajectories 
of X^(t) in (a) and of X^it) in (b). The fields behave stably; these reasonable results 
indicates that our numerical code works well. 

FIG. 3. Dynamics of a monopole for rj = 0.4mpi and A/e^ = 0.1. We assume two initial 
configurations: c$ = c„, = 1 in (a) and dp = Cw = 10 in (b). In (a) and (b) we plot the 
trajectories of X^{t) and X^(t) as well as apparent horizons. In (a) a monopole shrinks and 
the black-hole horizons appear, while in (b) a cosmological horizon exists from the beginning 
and a monopole continues to expand. We also draw the distributions of p + Spj in (c) and 
in (d), which correspond the results in (a) and in (b), respectively. In (c) the values around 
the center get negative at the beginning, but they bounce back to positive values, which 
confirms that the monopole core never infiates. On the other hand, in (d) the values of 
p + around the center remain negative from the beginning. This behavior indicates that 
exponential expansion really occurs inside the monopole. 
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FIG. 4. Dynamics of a monopole for = 0.55mpi, A/e^ = 0.1 and 0$ = Ci„ = 1. In (a) we 
plot the trajectories of X^{t) and Xy^{t) as well as apparent horizons. (Please also refer to 
Fig. 5, which is a schematic sketch of the spacetime structure.) From the beginning there 
are two apparent horizons, 5*1 and 5*2: 5*1 is the cosmological horizon. Later other two 
apparent horizons, S3 and 5*4, appear, and then 5*2 and 5*4 approach each other. These 
surfaces turn out to be black-hole horizons, 5*2' and 5*4', the moment they intersect. In (b) 
we draw the distributions of p + Sp^. The values around the center are initially positive, 
but they become negative. This suggests that, if t] is large enough, a monopole begins to 
expand exponentially even if its initial size is not so large. We also show in (c) the relation 
between the proper distance along the radial direction and the circumference radius, which 
indicates a wormhole structure really appears. 

FIG. 5. Schematic sketches of the spacetime structure. These figures are not generated 
from the numerical data, but they are based on the results presented in Fig. 4. 

FIG. 6. Dynamics of a monopole for r] = 0.3mpi and A/e^ = 1. Setting c$ = c^„ = 1, we 
show the evolution of $ in (a) and that of w in (b), and the trajectories of X^(t) and of X^^t) 
in (c). Although some oscillations remain outside the monopole, the core of the monopole 
approaches a stable configuration. We change the initial size in (d), finding monopoles with 
any initial size behave stably. These results indicate that there exist stationary solutions. 

FIG. 7. Formation of a monopole from the perturbed symmetric state (de Sitter spacetime). 
An initial configuration is given by ( ^.4| ). We draw the configuration of $ in (a) and that 
of w in (b), and the trajectories of Xq,{t) and Xy,{t) in (c). Once a monopole configuration 
is formed, the monopole approaches a stable configuration instead of continuing to expand. 
This result supports the existence of stationary solutions. 

FIG. 8. Parameter plane of r]/mpi and A/e^ in which we summarize our numerical results. 
A square (□) denotes a stable solution, as is the case in Fig. 2 or 6. A cross (x) denotes the 
case where a monopole shrinks, as is the case in Fig. 3(a). A circle (Q) denotes the case 
where a monopole inflates and the wormhole structure appears, as is the case in Fig. 3(c) or 
4. A dotted line indicates the maximum values of rj/mpi versus A/e^, depicted approximately 
by use of Fig. 6 in |Q. We vary c$ and Cw from 1 to 10, and hence some parameter points 
are labeled as two symbols. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9512045v3 



This figure "fig4c-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9512045v3 
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